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Randomized benchmarking of single qubit gates in a 2D array of neutral atom qubits 
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We characterize single qubit Clifford gate operations with randomized benchmarking in a 2D 
array of neutral atom qubits, and demonstrate global and site selected gates with high fidelity. 

An average fidelity of F 2 = 0.9983(14) is measured for global microwave driven gates applied to 
a 49 qubit array. Single site gates are implemented with a focused laser beam to Stark shift the 
microwaves into resonance at a selected site. At Stark selected single sites we observe F 2 = 0.9923(7) 
and an average spin flip crosstalk error at other sites of 0.002(9). 


PACS numbers: 03.67.-a, 31.10.Jk, 42.50.Dv, 03,67.Lx. 

Qubits encoded in hyperfine states of neutral atoms 
are one of several promising approaches for scalable im¬ 
plementations of quantum information processing[l|. In 
this letter we demonstrate and characterize single qubit 
gate operations in a 2D array of up to 49 neutral atom 
qubits encoded in long lived hyperfine states. Using a 
microwave field we implement arbitrary Bloch sphere ro¬ 
tations on either the entire 2D qubit array in parallel, 
or on single sites that are selected by an auxiliary Stark 
shifting laser beam[2|]. Single atom qubits are stochasti¬ 
cally loaded into the array, with an average of 29 sites 
filled for the data reported here. This is the largest num¬ 
ber of individually controllable qubits for which quantum 
gate operations have been characterized to date. The fi¬ 
delity of the global operations, site selected operations, 
and crosstalk during site selected operations are quanti¬ 
fied using randomized benchmarking (RB)i- 

Control of individual qubits in a spatially extended 
array is an important capability in ongoing efforts to de¬ 
velop scalable quantum processors. Atomic qubits en¬ 
coded in hyperfine ground states can be controlled with 
optical or microwave frequency fields. Optical fields 
can be tightly focused to address individual qubits as 
has been demonstrated in several experiments |4|. When 
the spatial separation of qubits is comparable to, or 
less than, the optical wavelength, addressing by focusing 
alone is not sufficient to suppress crosstalk to neighbor¬ 
ing sites. Addressing with subwavelength resolution can 
be achieved using quantum interference techniques [B| or 
by using an additional external field gradient to select a 
desired site. This latter method has been implemented 
with magnetic field gradients or with auxiliary Stark 
shifting optical beams in conjunction with microwave 
fields [2|, L7H9(. The use of microwave fields for qubit con¬ 
trol is particularly convenient since both the global rota¬ 
tions which are the starting point for many quantum al¬ 
gorithms, as well as single qubit control needed for gates, 
can be implemented with the same control hardware. 

In this letter we demonstrate global and site selected 
single qubit gates using microwave drive fields with a 
tightly focused Stark shift beam. We derive optimal val¬ 
ues for the Stark shift which minimize crosstalk to other 


sites. Arbitrary rotations on the Bloch sphere are im¬ 
plemented using variable length and phase microwave 
pulses. In contrast to previous experiments which used 
either adiabatic pulses that do not provide full control 
on the Bloch sphere[7j], or spatially periodic Stark shift¬ 
ing techniques which do not address single sites [§!, 9j, we 
demonstrate full control at single sites of a 2D qubit ar¬ 
ray. Using RB techniques we characterize the fidelity of 
Clifford groupjrntes, as well as the crosstalk during site 
selected gates[10[. 

RB was introduced in .3] as an efficient approach for 
characterization of quantum gate fidelities. It has sev¬ 
eral advantages compared to full tomography including 
resource requirements that scale linearly with the num¬ 
ber of qubits and the capability of distinguishing gate 
errors from state preparation and measurement (SPAM) 
errors. RB has been used for characterization of one- and 
two-qubit gates, as well as quantum processes, on a vari¬ 
ety of qubit platforms including ions[3l 11], nuclear mag¬ 
netic resonancejlii]. superconductors [l3| . neutral atoms 
[1B 3, and quantum dots|l6l- We encode qubits 
in the Cs clock states with |0) = \f = 3,to/ = 0), 
|1) = |/ = 4, to/ = 0). Our implementation of RB 
uses the complete set of 24 Clifford gates C\ . These 
are generated from the set {I, Rj/2), Rj(n)} where 
Rj (0) = e ~ l0crj / 2 with <jj Pauli matrices about axes 
j = x,y,z. We use constant amplitude pulses of mi¬ 
crowave radiation resonant with the uj q = 2 tt x 9.19 GHz 
|0) -SA |1) clock transition for R x rotations. Phase shifting 
the microwaves provides R y rotations. R z operations are 
implemented by composing x and y axis rotations. The 
microwave pulses used for each Clifford gate are listed in 
the supplementary material. 

The main elements of the atomic experiment are as de¬ 
scribed in 3. In brief, a two-dimensional array of blue 
detuned optical traps is defined using 780 nrn light pro¬ 
jected into a pyrex ultrahigh vacuum cell. The 7 x 7 = 49 
site array has a 3.8 p.m site to site spacing and trap 
depths of ~ 400 /iK for Cs atoms. The Cs atoms are 
collected in a 2D cooling region, transferred to the pyrex 
cell with a push beam, and then trapped and cooled in 
a 3D magneto-optical trap (MOT). The array is then 
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FIG. 1. (color online) Qubit array measurements, a) False- 
color fluorescence image of single atom qubits filtered with an 
independent component analysis. Each pixel views a region 
0.62 x 0.62 p m The image is a composite averaged over 500 
exposures, b) Histogram of the number of occupied sites for 
500 array loading measurements. The inset shows the single 
atom photoelectron counts at a single site for 2000 measure¬ 
ments. The average overlap of the Gaussian fits at all 49 sites 
was 0.0004. 


FIG. 2. (color online) a) Microwave Rabi oscillations at a 
single site during global addressing. Each point is the average 
of 50 measurements and the fitted Rabi frequency is H = 2-7T x 
4.74 kHz. b) Oscillations on site 31 using Stark addressing 
with the average probability of measuring a neighboring site 
in the wrong state due to crosstalk during a 37r oscillation 
indicated numerically. The Rabi frequency for this data was 
ft = 2n x 8.5 kHz and the sites are numbered from 0 in the 
upper left corner to 48 in the lower right corner. 


turned on, the MOT quadrupole field is turned off, and 
the captured atoms are cooled to 5-10 p iK using polar¬ 
ization gradient cooling. About 5 ms of near resonant 
light at 852 nm is used to invoke photo-assisted collisions 
which remove atoms in multiply occupied sites. After 
this preparation step multiple atoms are not observed at 
any site. The presence of an atom is detected by fluo¬ 
rescence imaging as shown in Fig[T| We observe single 
atom loading rates approaching 70% at a few sites, which 
is suggestive of repulsive light assisted collisions [18( . On 
average we load 60% of the array sites with a single atom 
in each experimental run. 

An initial fluorescence image reveals which sites are 
loaded with qubits for each experimental run. The atoms 
are then optically pumped into 11) using tt polarized 894 
nm light resonant with | 6 si/ 2 ,/ = 4) —> | 6 pi/ 2 ,/ = 4) 
and 7 r polarized 852 nm repump light resonant with 
| 6 si/ 2 ,/ = 3) —> | 6 p 3 / 2 ,/ = 4). The quantization axis 
is perpendicular to the plane of the array and is defined 
by a 0.15 mT magnetic bias field. After the quantum 
gate operations described below a state sensitive mea¬ 
surement is performed. To measure the probability of |0) 
we push out atoms in / = 4 with unbalanced resonant 
light pressure (| 6 si/ 2 ,/ = 4) —I | 6 p 3 / 2 ,/ = 5)) and then 
measure the presence of an atom by integrating the flu¬ 
orescence from MOT light with detuning —7^q P3/2 and 
resonant saturation parameter 3.3 for 20 ms. This re¬ 
sults in high fidelity discrimination of the qubit states as 
is seen in Fig. DO- To measure the probability of 11) 
we apply a R x ( 7 r) microwave pulse before the push out 
and fluorescence measurement. All sites in the array are 
measured in parallel using camera detection of the flu¬ 
orescence. The dominant error in state measurement is 
the small probability of transferring an atom from / = 4 
to / = 3 during the pushing out step. 

In preparation for quantum gate experiments the qubit 
array characteristics are measured. Array averaged val¬ 


ues are 17 s for the 1/e atom lifetime, 0.59 s for T) and 
14 ms for T 2 *. The T 2 value, which is measured using 
microwave Ramsey spectroscopy, is dominated by mag¬ 
netic noise and finite temperature motional effects 1 19|. 
The 9.19 GHz microwave source is locked to a GPS disci¬ 
plined crystal oscillator. The frequency is slightly shifted 
from the free space Cs clock transition due to the mag¬ 
netic bias field and the ~ 500 Hz light shift at the center 
of each trapping site. A maximum T 2 of ~ 50 ms has 
been observed at a few sites which we attribute to vari¬ 
ations in the cooling efficiency and atom temperature. 
We anticipate that the T 2 value can be substantially im¬ 
proved in future experiments using trap compensation 
techniques^. 

We proceed with RB experiments to measure the fi¬ 
delity of single qubit gate operations with microwave ra¬ 
diation from a horn external to the vacuum cell driving 
all qubits in parallel. Random Clifford gate sequences of 
length l are generated with each gate chosen uniformly 
from C\. The average pulse area per gate was 77 t /4 (see 
supplemental material). We start with all qubits in 1 1 ). 
At the end of each sequence we add a final gate which, 
in the absence of errors, should transfer the qubits to | 0 ). 
In the presence of depolarization errors the probability 
of measuring | 0 ) is 

+ (1) 

Here da is the depolarization probability associated with 
state preparation, the final transfer gate, and state mea¬ 
surement, while d is the average depolarization of a Clif¬ 
ford gate. Using the standard definition [2l[ of the fi¬ 
delity of two density matrices p , p' given by W(p, p') = 
Tr \\J\fpp'yfp\ one can readily show that the square of 
the average gate fidelity is F 2 = 1 — d/2. The quantity 
F 2 is equal to the average fidelity of a Clifford gate F a 
introduced in |3 . 
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FIG. 3. (color online) Probability of measuring the correct 
output state at site 27 of the array for 7 RB sequences. 
Each sequence was truncated at 10 different lengths l = 
{1,12, 23, 34, 45, 56, 67, 78, 89,100}. Each data point is an av¬ 
erage of 50 measurements. The inset shows a histogram of 
gate hdelities for 47 of the 49 array sites. Two sites were 
dropped due to poor loading statistics. 


We applied 7 randomized Clifford gate sequences for 
all the trapping sites. Representative data from a single 
site are shown in Fig. [3] together with a histogram of F 2 
across the array obtained by extracting d from fits to Eq. 
Jl]) at each site. The results are summarized in Table |U 
The highest fidelity seen at any site was 0.9999(3) with 
an array average of 0.9983(14). These values are compa¬ 
rable to the highest fidelity neutral atom gates reported 
to date[Hj], wh° reported a global average RB fidelity 
of 0.99986(1). An indication of where these experimen¬ 
tal results stand in relation to theoretical thresholds for 
fault tolerant quantum computing can be found by con¬ 
sulting Table 8 in Ref. 221. 

In order to understand the sources of the observed er¬ 
rors we simulated the RB experiment allowing for de¬ 
tuning from LOq by up to 100 Hz, which corresponds to 
the average differential Stark shift of the trapped qubit 
states across the array. We also included pulse length 
timing errors of up to 200 ns corresponding to 0.2% drifts 
of the microwave power. Accounting for these imperfec¬ 
tions predicts gate errors several times smaller than those 
observed. Including a density matrix coherence decay 
factor [li| a(t,T|) = l/2±(l/2)/[l±0.95(yr) 2 ] 3 / 2 we es¬ 
timate ( F 2 ) = 1— (d)/2 = 1 — [1— a((t)c 1 ,T£)\/2. Putting 
(t) Cl = (6)cJSl = (7tt/4)/(2t r x 4.74 kHz) = 185 ps 


and T 2 * = 2.7 ms we recover the observed (F 2 ) 4? sites = 
0.9983 from Tabled] The median T£ observed in the array 
is over twice longer at 7.0 ms. We conclude that the fac¬ 
tors limiting the gate fidelity found from RB experiments 
are a combination of finite T 2 * which could be improved 
using echo techniques or trap compensation]^, and er¬ 
rors in the pulse length and detuning. 

To perform gates on individual qubit sites we detuned 
the microwave frequency w by d = w — ~ 27r x 33 kHz. 

This detuning suppresses the microwave qubit rotation 


TABLE I. Results of RB fidelity measurements for global 
(first five rows) and single site addressing. The last three 
rows are (E xt ) the average crosstalk error on the entire array, 
(E x t)nn the average crosstalk for the nearest neighbor sites, 
and (E xt )_nn the array averaged crosstalk excluding the near¬ 
est neighbor sites. 


) 47 sites 

0.092 ± 0.066 

(d) 47 sites 

0.0035 ± 0.0027 

(F 2 ) 

\ / 47 sites 

0.9983 ± 0.0014 

F 2 in 

0.9939 ± 0.0007 

F 2 

max 

0.9999 ± 0.0003 

p 2 

-^single site 

0.9923 ± 0.0007 

(d[{ )xt 

0.037 ±0.027 

<£xt> 

0.002 ± 0.009 

(-E'xt)nn 

0.014 ±0.02 

(-Ext)-nn 

0.0005 ± 0.001 


by a factor scaling as f l 2 /5 2 . We then selected a de¬ 
sired site using a tightly focused 459 nm beam with 1 /e 2 
intensity radii of = 3.2, w y = 2.7 pm detuned by 
As = 27T x 20 GHz from the |6si/ 2 , / = 4) <->■ |7pi/ 2 , / = 
4) transition. The beam size was chosen as a compro¬ 
mise between tight focusing which gives small crosstalk 
to neighboring sites, and loose focusing which reduces 
sensitivity to beam misalignment on the target site. The 
intensity of the 459 nm beam was adjusted such that the 
induced differential Stark shift of states |0),|1) was set 
equal to 6 to bring a selected site into resonance Q. The 
459 nm Stark beam was cr+ polarized and propagated 
normal to the plane of the array. Using a pair of orthog¬ 
onal acousto-optic modulators the Stark beam could be 
scanned to a desired qubit site with a switching speed 
under 0.5 ps. 

The choice of detuning <5 for single site gates is a 
trade off between less than perfect suppression of the mi¬ 
crowave field at small <5 and excessive photon scattering 
from the Stark beam at large 6. For a given value of <5 
the photon scattering can be reduced by working at large 
optical detuning As, but not completely eliminated since 
for large As, 8 ~ ui q / A| and tends to zero. An optimized 
working point which reduces the need for large S can be 
found by choosing a detuning for which the off-resonant 
coupling to non-selected sites gives a pulse area which is 
a multiple of 4-7 t and therefore does not disturb the qubit 
states. For a pulse area of Or on the targeted qubit the 
condition for minimal disturbance of non-targeted sites 
is 8/fl = {n 2 \d)TT 2 / 6 2 r — l) 1 ^ 2 with n an integer. Thus 
the leakage error should have a first local minimum for a 
7r pulse at (5/fl ~ \/l5. 

This estimate can be verified by a calculation which av¬ 
erages over all possible states of the non-targeted qubits. 
Let the initial state be | Q,<j>) = cos(|)|0) + e 1 ^ sin(|)|l). 
This state receives a unitary transformation \6, <j>) — > 
Uj{6 r, 5)\9, (j)) with Uj{6 R ,8) the operator for a 9r ro- 
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FIG. 4. (color online) Crosstalk error for JU(7r/2)(blue solid 
line), Fix (rr)(yellow dashed line), _R z (-7r/2)(blue dotted line), 
and R z (7r) (yellow dashed-dotted line) rotations. The vertical 
dashed lines mark values of y/16n 2 — 1. Panel a) is the error 
from Eq. (2) averaged over input states and panel b) is the 
error for the initial state 11). 


tation about axis j detuned by S. The fidelity of the 
transformed state with respect to the original state, av¬ 
eraged over the Bloch sphere, is 

^ _ fidO sin(g)/ o ^#|(0,^|f/ J (^,J)|g,^)| 2 

^ 4 7T 

( 2 ) 

The crosstalk error defined as E xt = 1 — F(6 r, 6) is shown 
in Fig. 2] for several elements of C\. We see that the 
simple estimate of <5/f2 ~ \Zl5 for a R x {n) rotation is 
verified by the full calculation. Since the error is min¬ 
imized at different detunings for different rotations the 
detuning should be dynamically adjusted in tact with 
the gate being performed. In the demonstration de¬ 
scribed below we have simply used a fixed detuning of 
^= 2 2 ;x X 8 33 5 k k H Z =3.88cVl5. 

To characterize site selected gates we applied 10 ran¬ 
domized Clifford sequences as shown in Fig. [5] Each 
sequence was truncated at 8 different lengths {1, 8, 15, 
22, 29, 36, 43, 50}. Averaging over the 10 sequences 
yields an average gate fidelity F 2 = 0.9923, giving a per 
gate error which is about 4.5 times larger than for the 
array averaged global gates. We attribute this to fluctu¬ 
ations in the intensity and pointing stability of the Stark 
shifting beam resulting in deviations from the optimal 
detuning condition. 

The crosstalk error at other sites was measured by 
preparing them in |1) and then measuring P| 0 ) = ^ — 



number of gates 

FIG. 5. (color online) RB data at Stark beam addressed site 
31 for 10 RB sequences. Each data point is an average of 50 
measurements. 


^(1 — dif)(l — d x tY after each Clifford sequence was ap¬ 
plied. Dropping the Stark addressed site and sites whose 
loading was poor yields an average (d x t). The array aver¬ 
aged background error on non-addressed sites per Clifford 
was (E xt ) = (d x t) /2 = 0.002(9). Due to the finite size 
of the Stark beam there was intensity overlap to near¬ 
est neighbor sites that was as high as 5%, resulting in 
increased crosstalk compared to further away sites. The 
crosstalk values for the neighboring sites are given in Ta¬ 
ble m The average crosstalk error is comparable to the 
average error of global gates. However, this result was 
obtained for an initial state in the non targeted sites of |0) 
and therefore essentially corresponds to a spin flip error. 
It is to be expected that measurements with arbitrary ini¬ 
tial states would yield higher errors (compare Figs. |4ji) 
and b)). Ultimately, a slightly smaller Stark beam, and 
larger values of (5/ST than have been demonstrated here 
should be used for effective crosstalk suppression. 

In summary we have demonstrated high fidelity sin¬ 
gle qubit gate operations in a 2D array of neutral atom 
qubits. Using microwave pulses we perform either paral¬ 
lel gates on all qubits, or gates on single qubits selected 
by a Stark shifting beam. The results reported, together 
with the demonstration of two-qubit entanglement in the 
array using Rydberg blockade gates, which we will report 
on elsewhere [23], are a step towards scalable quantum 
computing with neutral atom qubits. 

This work was supported by the IARPA MQCO pro¬ 
gram through ARO contract W911NF-10-1-0347. 
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SUPPLEMENTARY MATERIAL FOR RANDOMIZED BENCHMARKING OF SINGLE QUBIT GATES 

IN A 2D ARRAY OF NEUTRAL ATOM QUBITS 


In this supplementary material we specify the microwave pulses used for each of the 24 gates in the single qubit 
Clifford group G\. The group elements and our implementation are listed in Table [TJ Note that rotations of —7r/2 were 
implemented as +37r/2 rotations. We therefore used microwave pulses which were longer than necessary giving an 
average pulse area per Clifford gate of (0)ci = 77 t/4 . If all 37 t/2 operations had been replaced by —7r/2 rotations this 
would have reduced the average pulse area to 557 t/48. As discussed in the main text the finite T% time accounts for 
about half of the observed gate error. This could be improved on simply by using the shorter —7 t/2 rotations instead 
of 37 t/2 rotations. 
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TABLE I. Elements of the Clifford group Ci for a single qubit. The notation 9 in column 2,3 or 4 is shorthand for Rj(9) about 
the corresponding axis. The resulting operator U is due to the sequence U X U V U Z , reading right to left. The operators have 
been factored so that the first nonzero element in the top row is unity. The microwave implementation R 3 R 2 R 1 , reading right 
to left, used the rotations shown in the last three columns. These give the same Clifford operations as U, up to irrelevant global 
phase factors. 





















